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We investigate the inhomogeneous chiral dynamics of the 0(4) linear sigma model in 1+1 dimen- 
sions using the time dependent variational approach in the space spanned by the squeezed states. 
We compare two cases, with and without the Gaussian approximation for the Green's functions. We 
show that mode-mode correlation plays a decisive role in the out-of-equilibrium quantum dynamics 
of domain formation and squeezing of states. 



The possibility of the formation of the disoriented chi- 
ral condensate (DCC) in high energy heavy ion collisions 
has been extensively studied with various methods. In 
classical approximation 0, , it has been shown that the 
amplification of long wavelength modes of the pion fields 
takes place when the system starts with the nonequilib- 
rium initial condition, quench initial condition In 
addition to the amplification, spatial correlation of the 
fields has been also shown to grow. 

Although the classical approximation is expected to 
work well in incorporating nonequilibrium aspects of the 
system when pion density is large, it is still desirable to 
include quantum effects. In fact, investigations in this 
direction have been also carried out extensively with the 
Hartree approximation, the large N approximation, and 
so on 0,01 . In most of the previous studies which include 
quantum effects, however, it has been assumed that the 
system is spatially homogeneous. Problems such as in- 
sufficient thcrmalization at late times and impossibility 
to describe domain structures have been recognized. It 
has not been conclusive whether there is a chance for the 
correlations to grow through nonequilibrium time evolu- 
tion. 

There are at least two ways for possible improvement. 
One is to include higher order quantum corrections and 
the other is to accommodate spatial inhomogeneity. We 
will pursue the latter in this paper. Recently, the dynam- 
ics of spatially inhomogeneous system has been studied 
by several groups quantum mechanically [MS 00 an d it 
has been shown that the thermalization of the quantum 
fields can occur. In these works, the Gaussian approxi- 
mation, in which the Green's functions are assumed to be 
diagonal in momentum space, has been adopted because 
of computational reasons. Physically, it corresponds to 
ignoring correlations between modes with different mo- 
menta, and under the approximation different modes can 
interact only through the mean fields. However, it is pos- 
sible that the direct coupling of modes through the off- 
diagonal correlations is important for the time evolution 
of the system when the system does not possess transla- 
tional invariance. To see if such an effect is substantial, 
we study the dynamics of chiral phase transition in spa- 
tially inhomogeneous systems with off-diagonal compo- 



nents of the Green's function in momentum space fully 
taken into account. 

In this paper, we take the 0(4) linear sigma model as a 
low energy effective theory of QCD and apply the method 
of the time dependent variational approach (TDVA) with 
squeezed states. This method was originally developed 
by Jackiw and Kerman as an approximation in the func- 
tional Schrodinger approach |9| and later it was shown 
to be equivalent to TDVA with squeezed states by Tsue 
and Fujiwara [Tof. 

In this approach, the trial state is a squeezed state 
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where a runs from to 3. a = is for the sigma field and 
a = 1 — 3 are for the pion fields. |0) is the reference vac- 
uum and G^\x,y) = (O|0 Q (f)^> a (y)|O). (f> a {x) and ir a (x) 
are the field operator and conjugate field operator for the 
field a, respectively. C a (x, t), D a (x,t), G a (x,y,t), and 
II a (x, y, t) are the mean field variables for (j) a field at x 
and t, the canonical conjugate variable for the mean field, 
the quantum correlation for x ^ y (fluctuation around 
the mean field for x — y), and the canonical conjugate 
variable for G a (x, y, t) , respectively, and all of them are 
real functions. N a (t) is a normalization constant. T(t) is 
an operator that describes the squeezing, and if T a (t) is 
set to 0, |<& a (i)) is reduced to a coherent state with the 
expectation values of 4>(x, t) and tt(x, t) given by C a (x, t) 
and D a (x,t), respectively. 

The Hamiltonian H of the 0(4) linear sigma model is 
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given by 
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As shown in Eq. (JTJ, the trial state is specified by 
C a (x,t), D a (x,t), G a {x,y,t), and R a (x,y,t). Their time 
evolution is determined through the time dependent vari- 
ational principle: 



-,Jdtm)\i§£-Hm))=o- 



(3) 



In this approach, correlation between different modes in 
momentum space arises through the scattering of quanta 
caused by the nonlinear coupling term in the model 
Hamiltonian even if there is initially no such correlation. 
This can be seen best from the following equations of 
motion in momentum space, 
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where m 2 = 4Au 2 , and C a (k,t), G a (k,k',t), and 
H a (k,k' ,t) are the mean fields for the </> a field with mo- 
mentum k, the correlation between modes with momenta 
k and k' for k ^ k' (the quantum fluctuation around the 
mean field for k = k'), the canonical conjugate variable 
for G a (k,k' ,t), respectively. In Eq. Q, we have used 
the notation, 
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dk" 
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In the Gaussian approximation, G a (k,k' \t) and 
n a (fc, k', t) are set to zero for k ^ k' and correlations be- 
tween different modes in momentum space are ignored. 
However, MS 2 '(h — k',t) in Eq. (0J, which originates 
from the four-point interaction terms in the 0(4) linear 
sigma model, couples modes with different momenta and 
correlations between them develop even if initially there 
exists no correlation among them. 



In numerical calculation, we have assumed the one- 
dimensional spatial dependence for the mean fields and 
the Green's functions for computational simplicity. In 
addition, we have imposed the periodic boundary condi- 
tion for the mean fields and the Green's functions. We 
have carried out calculation on a lattice with the lattice 
spacing d = 1.0 fm and the total length L = 64 fm, which 
leads to the momentum cutoff A = 1071 MeV. The pa- 
rameters A, v, and ho are determined so that they give 
the pion mass M n = 138 MeV, the sigma meson mass 
M„ = 500 MeV, and the pion decay constant = 93 
MeV in the vacuum following the prescription given in 
Ref. [13, and we have obtained A = 3.44, v = 110 MeV, 
and h = (103 MeV) 3 . 

There are several scenarios for the DCC formation in 
high energy heavy ion collisions. Here we adopt the 
quench scenario. In this scenario, the chiral order pa- 
rameters remain around the top of the Mexican hat po- 
tential after the rapid change of the effective potential 
from the chirally symmetric phase to the chirally broken 
phase. In order to take this situation into account, we 
have used the following initial condition; at each lattice 
point, the mean field variable for the chiral fields and 
their conjugate variables C a (x, 0) and D a (x, 0) are ran- 
domly distributed according to the Gaussian form with 
the following parameters [13] , 
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where D = 1 is the spatial dimension and 5 is the Gaus- 
sian width. We shall use S = 0.19v in the following calcu- 
lations. In relating the Gaussian widths of C a (x, 0) and 
D a (x, 0) , we have taken advantage of the virial theorem 
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As for the initial conditions for the quantum fluctua- 
tion and correlation, we have assumed that their values 
are those realized in the case where each state in mo- 
mentum space were in a coherent state with a degenerate 
mass toq for the sigma meson and pions, namely, 



G a (x,y,0) = 



n a (f,y,0) = 0, 



dk 1 
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where LUk — yml + k 2 . We adopt to = 200 MeV. 
As shown above, the Green's functions G a (a;, y, 0) and 
IL a (x, y, 0) are initially diagonal in momentum space. 
The off-diagonal components appear in the course of the 
time evolution of the system due to the direct mode-mode 
coupling induced by A^ 2 ) in Eq. (@J. 

We have carried out two sets of numerical calculations. 
In one case, we have taken into account all components 
of the two-point Green's functions (case I), while in the 
other case only the diagonal components of the two-point 
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FIG. 1: (Color online) Time evolution of the mean fields 
in case I. (a) and (b) are for the sigma field and the third 
component of the pion field, respectively. 



FIG. 2: (Color online) Time evolution of mean fields in case 
II. (a) and (b) are for the sigma field and the third component 
of the pion field, respectively. 



Green's functions (case II) were included in the calcula- 
tion as in most of the preceding works. 

In Fig.^ we show the time evolution of the mean fields 
for the sigma and the third component of the pion fields 
(</>o and </>3, respectively) obtained with the calculation 
including all components of the Green's functions (case 
I). As a whole, the expectation value of the sigma field ap- 
proaches a constant. On the other hand, that of the pion 
field oscillates around zero and shows a domain struc- 
ture with long correlation length. This is the formation 
of DCC domains. It is observed that the domain struc- 
ture continues to grow till as late as t ~ 40 fm. 

In Fig. [21 we show the time evolution of the same 
mean fields obtained with the Gaussian approximation, 
i.e., without the off-diagonal components of the Green's 
functions (case II). At the beginning, the behavior of the 
mean fields is similar to that in case I. However, after 
a few femtometers, there appears a clear difference be- 
tween the two cases. In case II, short range fluctuation 
is dominant and no long length correlation is observed. 
No qualitative change in the behavior of the mean fields 
takes place in case II after a few femtometers. This tells 
us that the mode-mode correlation plays a decisive role 
in the formation of DCC domains. 

To examine the growth of spatial correlation of the 
pion fields more quantitatively, we define the following 



spatial correlation function C(r, i), 

c , r a = I $(x)-<f(y)S(\x-y\ -r)dxdy 
I \^>{x)\\(l>{y)\S(\x - y\ ~ r)dxdy 

w here $(x) ■ $(y) = £\ =1 4>i(x)(j>i{y) and \<j)(x)\ = 

In Figs. Ola) and[3fb), we show this spatial correlation 
function in cases I and II, respectively. The correlation 
functions are calculated by taking average over 10 events. 
Substantial generation of the correlation takes place af- 
ter the typical time scale of the initial rolling down of the 
chiral fields, say, a few femtometers in case I, while the 
growth of the spatial correlation ends in case II by t ~ 5 
fm. The domain formation of DCC shown in Fig. ^b) 
and this growth of spatial correlation shown in Fig. EJa) 
beyond the rolling down time scale may be related to the 
parametric resonance Q ■ We are currently investi- 
gating this possibility. 

Next we show the time evolution of the quantum fluc- 
tuation, which is represented by the same-point Green's 
function G a (x,x,t). We define the spatially averaged 
fluctuation function (F a (t)) spa cc at time t as 

(F a (t)) spax;e = ^ J G a (x,x,t)dx, (9) 
where V is the volume of the system. 
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FIG. 3: (Color online) The spatial correlations for the pion 
fields in case I (a) and case II (b). 
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FIG. 4: (Color online) The time evolution of the spatially 
averaged quantum fluctuation of the third component of the 
pion field, {Fs(t)) apa , ce for cases I and II. 



In Fig. we compare the time evolution of the spa- 
tially averaged quantum fluctuation of the the third com- 
ponent of the pion held, (i 7 3(i)} S pacc in the two cases. 
We observe remarkable increase of quantum fluctuation 



in case I, while only small amplification is seen in case 
II. This increase also lasts until about t = 40 fm. The 
comparison between case I and case II tells us that the 
off-diagonal correlations take an important role also for 
the enhancement of the quantum fluctuation. We have in 
fact found the growth of off-diagonal components of the 
Green's functions in momentum space. Note that this 
phenomenon cannot be described unless the off-diagonal 
correlation is introduced. This has an important meaning 
also for the identical particle correlation in high energy 
heavy ion collisions. Usually it is assumed that identical 
particles with different momenta are emitted indepen- 
dently. However, if there is quantum correlation between 
two different modes, this assumption becomes invalid, 
and it will be necessary to reformulate the theory for the 
identical particle correlation. 

In summary, we have studied the inhomogeneous chi- 
ral dynamics of the 0(4) linear sigma model in 1+1 di- 
mensions using TDVA with squeezed states. We have 
compared two cases. One is a general case in which both 
the mean fields and the Green's functions are inhomo- 
geneous, and the other is a case with the Gaussian ap- 
proximation, where translational invariance is imposed 
on the Green's functions. We have shown for the first 
time that the large correlated domains can be realized in 
the quench scenario with quantum mechanical treatment. 
More specifically, we have shown that the large amplifica- 
tion of quantum fluctuation and large domain structure 
emerge when all components of the Green's functions are 
retained, while only small quantum fluctuation and small 
and noisy domain structure are seen in the case with the 
Gaussian approximation. The dynamics in 1 + 2 and 1 + 3 
dimensional cases is of great interest and indispensable 
for the understanding of the DCC formation in ultrarel- 
ativistic heavy ion collisions. We expect more enhanced 
domain growth in 1 + 2 and 1 + 3 dimensional cases be- 
cause of the existence of more mode-mode correlations 
in such cases. We plan to confirm it by actual numerical 
calculation. 
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